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TECHNICAL NOTE 2350 , 


ON THE SECOND-ORUER TUNNEL-WALL-CONSTRICTION -CORRECTIONS 
IN TWO-DIMENSIONAL COMPRESSIBLE FLOW 
By E. B. Klunker and Keith C. Harder 


SUMMARY 


Solutions of the first- and second-order Prandtl -Busemann iteration 
equations are obtained for the flow past thin, sharp-nose, symmetric, 
two-dimensional bodies in closed channels. With the use of these solu- 
tions an expression is derived for the tunnel- wall interference. The 
tunnel-wall correction for a parabolic-arc airfoil is calculated to 
indicate the effects of compressibility, ratio of the tunnel height to 
the airfoil chord, and airfoil thickness coefficient. It appears that, 
for cases where the tunnel-wall corrections are significant, both the 
second- order effects and the variation of the correction along the chord 
should be considered. 


INTRODUCTION ' 


The use of wind tunnels for solving many aerodynamic problems in 
the high-subsonic speed range makes it desirable to reexamine the ques- 
tion of tunnel-wall corrections. The general problem of tunnel-wall 
interference in inc oirpre s s ible- f low fields has been treated in refer- 
ence 1. There the various types of interference are described and 
equations are presented which permit the correction of wind-tunnel data 
to free- flight conditions. These results may be carried over to include 
the first-order effects of compressibility with the aid of the familiar 
Prandtl-Glauert rule. (See reference 2, for example. ) It is well-known, 
however, that the linearized compressible-flow relations in general do 
not describe the phenomena accurately at high-subsonic Mach numbers. 

The exact analytical solution for the compressible flow past an 
arbitrary closed body has not been formulated; however, certain approxi- 
mate methods prove useful for studying flow problems. The Rayleigh- 
Janzen method, in which the initial step is the complete incompressible- 
flow solution, has been used frequently for obtaining approximate solu- 
tions to aerodynamic problems. In reference 3 the flow past a circular 
cylinder in a tunnel has been obtained by this procedure. The Rayleigh- 
Janzen method is -best suited, however, for obtaining flows past relatively 
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“hi nnt. profiles at low Mach numbers; hence, this method is not suitable 
for evaluating the tunnel-wall interference for thin profiles at high- 
subsonic Mach numbers. Solutions of flow problems by relaxation methods 
have proved useful in many instances; this approach has been employed 
in references 4 and 5 to evaluate the tunnel- waH-constriction correc- 
tions for the NACA 0012 airfoil and the Kaplan bump, respectively. 
Although accurate solutions can be obtained by this method, lengthy 
numerical calculations are involved and the result is confined to the 
specific profile and channel configuration treated. 

The Prandtl -Euse mann small-disturbance iteration method has been 
found most useful in the evaluation of the flow over thin bodies in an 
infinite stream. This method therefore is employed herein to evaluate 
the second-order t unn el-wall corrections for the flow past thin sym- 
metric bodies at zero incidence. The third- and higher-order solutions 
can also be obtained by this method although the labor involved becomes 
prohibitive. 

The wall-interference effects are often calculated from a suitable 
system of images. Other approaches are possible, however, and in many 
instances prove more convenient. In particular, an integral representa- 
tion proves useful for evaluating some interference effects. This repre- 
sentation is employed herein to evaluate the wall interference (con- 
striction effects) at subsonic speeds for a thin, two-dimensional, sym- 
metric, sharp-nose body in a tunnel. The method is closely related to 
the familiar source-sink concepts for evaluating the flow over a thin 
body in an unbounded stream. The extension of the source-sink concepts 
to the solution of the second-order Prandtl-Busemann iteration equations 
has been previously discussed in reference 6. The method of solution of 
the interference problem to the second order thus parallels this work. 

As an example of the use of the equations developed herein, the wall- 
interference effects are evaluated for a symmetric parabolic-arc airfoil 
in a two-dimensional channel. 


ANALYSIS 


The exact nonlinear equation governing the two-dimensional flow of 
a compressible fluid for the velocity potential is J 
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■where 

x’, y* rectangular Cartesian coordinates in the flow plane 

u’, v* velocity components along the x’- and y*-axes, respectively 

a local speed of sound 


and the subscripts x* and y* denote partial differentiation with 
respect to these designated variables. With the introduction of a 
characteristic length c/ 2 , where c is' the chord, and the undisturbed 
stream velocity U as the unit of velocity, the potential equation may 
be written in the nondimens ional form 


C 1 




3>. 


xy 


£ 



(1) 


where x 





D = 


23> t 

Uc 


The Prandtl-Busemann iteration equations are developed from the 
nonlinear potential equation (l) (see reference 7 , for example) by 
assuming that <& may be expanded in the form 3> = x + + • • • 

where 0 n+ p and its derivatives are small compared to 0 n and 

is of the oider t n where t is the thickness coefficient of the air- 
foil. With these assumptions the first two equations of the Prandtl- 
Busemann iteration method are 


P^lxx + ^lyy = 0 (2a) 




where a 


y + 1 0 p 

= — — k-, B = 1 - Moo , is the stream Mach number, and 


2 p2 


7 is the ratio of specific heats at constant pressure and constant 
volume. 
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The constraint of the tunnel walls, or the tunnel-wall correction, 
is expressed as the difference between the velocity at any point in the 
flow in the tunnel and the velocity at the same point in an unbounded 
stream. For aeronautical applications this velocity increment is of 
interest primarily on the surface of the body, that is, for correcting 
the surface pressures on a body in a wind tunnel to free-stream condi- 
tions. This tunnel correction on the surface of the body is found to 
the second order by solving equations (2) subject to proper boundary 
conditions . 

Throughout this paper a bar over a quantity (for exanple, G>) 
represents the quantity in the channel whereas the same quantity with- 
out the bar (for example, $) denotes the quantity in an unbounded 
stream. The differential equations (2) are, of course, the same for 
the flow in the channel and in a free stream. ' 


Boundary Conditions 


Let the equations 


Y = tY 1 (x) 


and 


Y = 0 


(-1 < x < 1) (3a) 

( X £:L; X <-1) (3b) 


define a thin symmetric body of thickness coefficient t lying on the 
x-axis between -1 and 1. Then the boundary conditions for the flow 
over the body in an unbounded stream are: 

At infinity 


$x ~ ^ } ®y ■ ® 


(4a) 


and on the body 


$ y (x,Y) = Y’3 > x (x,Y) 


(4b) 


where Y* denotes the slope dY/dx. For the flow over the body defined 
by equations (3) in a tunnel (fig. l) whose walls are at y* = ±h* (or 
y = ±h) the boundary conditions are: 
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At the tunnel wall 

(5a) 
(5b) 

and at x— 

$ x = 1 , = 0 (5c) 

The boundary conditions for 0 -j_ and 0 2 to the order t and t 2 , 
respectively, are: 

At infinity 

01x = 01y = 0 (6a) 

on the body 

0 ly (x, 0) « tY» (6b) 

1 

and at infinity 


®y(x,h) = 0 

on the body 

$ y (x,Y) = Y»® x (x,Y) 


02x “ 02y " 0 (Ta) 

on the body 

02y(x,O) = tY 1 «0 lx (x,O) + P 2 tY 1 0 lxx (x,O) (7b) 

Similarly, the boundary conditions for 0^ and 0 2 to the order t 
and t 2 , respectively, are: . 

As X — ± co 


01x “ 01y = 0 


(8a) 
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at the wall y = h 

0 ly (x,h) = 0 (8b) 

on the body 

0 ly (x,o) = tYi* (8c) 

and as x— > 1“ 

02x = 02y = 0 (9a) 

at the wall y = h 

0 2y (x J h) = 0 ( 91 j) 

on the body 

1 

02y( x J°) “ tY l*01x( x >°) + P^i&jcxfoO) (9c) 

The solution of the nohhomogeneous second-order equation (2b) may 
be expressed as the sum of a particular integral ^(xjY) (or ^(^y)) 

and a function y) (° r ?2( x >y)) satisfying the homogeneous 

equation 


P ?( P2xx + c P2yy = 0 

The particular integral for equation (2b) (reference 8) is 


*2( x >y) = ^^lxj^ 1 + 5)01 " f Y01y] 


( 10 ) 


(ID 
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Then from equations (7) and. ( 9 ) the boundary conditions for cp 2 and. 

qvj are ‘ 

At infinity 

<P2x = c P2y = 0 

on the body 

q>2y(x,°) = F(x) 

where 

F(x) = P^tY^’02_x(x, 0) + P^tY~| yy (x, 0) - 4 - ^^t0i(x, O)0] J{ y(x,O) 

and. at x — >±00 

~^ 2 x = ? 2 y = 0 (12a) 

at the wall y = h 

? 2 y(^h) = f(x) (12b) 

on the body 

q>2y(x,0) =F(x) (12c) 

where 

F(x) = P 2 tY 1 «0 lx (x,O) + P 2 tY 1 0 1v ^(x,O) - Mj^l + ^)t0^(x,O) j^yUjO) 


f(x) b M^ 2 g h0 1x (x,h)0 lYY (x,h) 
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The tunnel-wall interference to the second order then is found from the 
solution of equations (2) subject to the, given boundary conditions. 


Solution of the Boundary-Value Problem 

First-order solution .- A solution of the first of equations (2) 
for 0^, satisfying the first two boundary conditions of equations (8), 
is, for y ^ 0, 


0l(x,y) 



cosh Pa(y - h) 
sinh pah 


cos a(x - I ) d£ 


from which 


_ poo poo 

01y( x '°) = Jcj. ^ J g( i)cos a(x - £) d£ 


(13) 


The, body slope tY^’ may he represented by its Fourier integral as 

d£ (14) 


tYi’(x) = — f da f Yn*(i)coB a(x - i) 
n J 0 


From equations (13) and (l^) the boundary condition for tangential flow 
(equation (8c)) on the body is satisfied by choosing g(l) = ~ Y^* (S). 
Thus the solution for 0q(x,y) is 


M*,y) - - h f v / 1 *i’(» b) cos a(x - 6 ) 

x P^o a sinh Pah 


d£ 


Inverting the order of integration and integrating with respect to a 
gives 


0l(x,y) = f Y]_* (i)loge 


2X 


2 jcmsh 2X(x - I) - cos 2ipy 


r & (15) 
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■where X, = -2— . 

2ph 


The solution 0 1 (x,y) satisfying the "boundary condi- 


tions given by equations (6) may be found from equation (15) by letting 
h approach infinity. Thus 


0!(x,y) “ lim 0 1 (x,y)=^ ? r Y x ‘ ( | ) l 0ge ^(7317^7^2 dl (l6) 

h — >00 P' 1 » 


■which is the well-known solution for the flow past a symmetric body in 
an unbounded stream by a distribution of sources. 

With Anq = 01 x( x jO) “ 01 v(x, 0 ) , the difference in velocity on the 
body to the first order is 

^1 = Yl '^^ coth ^( x - &) - - - 1 ^ jdl ( 17 ) 

Th£ first-order tunnel-wall correction to the x-component of 
velocity on the body is given by equation (17). This solution was 
obtained by considering the flow in the upper half-plane; the solution 
in the lower half-plane is known by symmetry. 

Second-order solution .- The second-order solution for the flow over 
a symmetric body by an extension of the concepts of source-sink distribu- 
tions has been discussed in reference 6. The second-order solution 0 2 
is given by 02 = 92 + ^2 "wh er e 92 satisfies the homogeneous equa- 
tion (10) and i[r 2 is the particular integral given by equation (ll). 

The form of the differential equation and boundary conditions for cp 2 

are the same as for 0 j_. Thus the solution for 0 2 may be written 
immediately as 

02 = F(|)loge ]j(x - I) 2 + d| + ijr 2 

and 

02 x( x >°) + ^(l + f)[ 0 l( x ^O) 0 ixx( x ^ 0 ) + 01x 2 (' x ,o3 

(18) 
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(It should "be noted that the second-order solutions given here are not 
valid for blunt-hose profiles since for these cases the particular 
integral introduces singularities /which cannot he canceled by a source 
distribution. ) 

The second-order solution for the flow over a symmetric body in a 
channel may be found in a manner analogous to that for Consider 

the expression 



cosh Ba(y - h) 
sinh Pah 


cos a.(x - I) d£ + 


1 

P* 



f(l) 


cosh pay / „ \ , t 

— 1 cos a(x - |) d£ 

sinh Pah 


(19) 


which is a solution of the homogeneous differential equation (10). Since, 
from equation (19) , 


<P2y( x >°) 




cos a(x - I ) dl 


and 


cp^x^) 




f(£)cos a(x - |) d| 


which equal F(x) and f(x), respectively, by the Fourier integral 
theorem, the boundary conditions of equations (12) on the body and at 
the wall are satisfied. Inverting the order of integration of equa- 
tion (19) and integrating with respect to a gives 


cp 2 


F(|)loge 


5F& 


= — f 

2p* l 

from which the boundary conditions at 


2l(x - 0 - cos 

2X(x - I) + cos 
infinity may be 



shown to be satisfied. 
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The quantity 02x( x ->^) ^ oun ^ after some reduction as 

02x( x j°) = p7 F U) coth X ( x - 6) d£ - 

^ ^ fU) tanh X(x - £) d£ + t2x( x >°) 
* ^ —00 


(20) 


The second.-ord.er correction 'to the Telocity on the surface of a 
body in a wind, tunnel is 

Au = Auq + Au 2 + 0(t 3 ) 


where 


A&i = 0ix( x ^°) " 01x( x j°) 


Au 2 = 02x( x >°) ” 02x( x >°) 


or from equations (17), (l8), and (20) 


Au 


- i hf_ x y i ,( i) £ <=° tt ^ + 

— f F(|) coth X(x - £) d£ - ~ f f(£) tanh \(z 
P*J_1 P«0_ w 


- I) d£ - 


r(5) -^ r + 

X - 5 


M^l + f)j0l(x,O)0 lxx (x ; O) +7lx 2 (x,0) - 0i(x,O)0ixx(x,O) 

fe’W,] + 0(t 3 ) 


( 21 ) 
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Throughout this paper the Cauchy principal value of all improper integrals 
is to he taken. 


Discussion of the Equations 


The tunnel-wall- interference effects are often evaluated hy con- 
sidering a suitable system of images. (See references 1 and 2, for 
example. ) The first-order equations developed previously are consistent 
with this approach. An expansion for coth X(x - £ ) (see reference 9 > 
for exanqple) is 


coth X(x - |) 


X(x - |) 


+ 2 


V” *•(* ~ 6) 

n=l n 2 it 2 + X 2 (x - 


£) 2 


Then, from equation (17) , 


2Xt 

AUl = piT" 



*•(* - I) 

X 2 (x - |) 2 + n 2 it 2 


d£ 


which is recognized as the influence of source-sink images located at 
y = The series appearing in the expression for Au^ is uniformly 

convergent; hence the term-by-term integration is valid. 


The calculation of the tunnel-wall interference for a given profile 
requires the evaluation of several integrals. The integrals occurring 
in equation (21) are not easily evaluated for most profiles and it is 
preferable to expand the integrands in the form of a series. The series 
given previously, however, is not the most useful. The expansions 


-,2n-v 2 u t 


X coth X(x - £) = — ^ — + / (-l) n (x - |) 

x - I n =l (2a) J 


2n-l 


tanh 


(X2(x - &)2 < * 2 ) 

X(x - I) = - 1)Xai ' lB n ( X . 


n=l 


(2n)J 


(x 2 (x - I ) 2 < £) 


( 22 ) 
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(see reference 10, for example), where the coefficients Bq are the 
Bernoulli numbers, are better suited for calculation since the integra- 
tions are more easily performed. Moreover these alternating series 
converge rapidly and in most cases only a few terms need be retained. 


With the use of the relations (22) the first-order correction to 
the velocity on the body is determined from equation ( 17 ) as 


Aui 


1 = ^ 2Z (-D 


n=l 


n-1 2 2a X 2n B n 

(2a) I 


-1 


V(6)(* " ^) 2n " 1 ^ 


(X 2 (x - |) 2 < it 2 ) ( 23 ) 


or, after integration by parts. 


4u x - 1 y ~ (-D n_i ( _ ~ - i)b p r 1 y 1 (d(x . s)2“- 2 m 

9 * (2a)! J -1 


U 2 (x - I) 2 < * 2 ) 

Thus the first term is proportional to the area of the body and the 
remaining terms are in a form convenient for graphical integration. 

Similarly, to the second order, the correction to the surface 
velocity is 


Au 


sr*-ht + 

«=i (2a) l 


n=l 


i f- (.if 1 r 1 ?(!)<x - + 1 r mi - m dl 

u* -Sr (2n) ; J.1 PxJ.l x-i 


i 2: (-D 0 - 1 g£n(gai - ^ r ♦ 

P* hr (2n): J -1 


l2Ut 


Mco 2 (l + f)|0i(x,O)^ lxx (x,O) + 0ix 2 U,O) - 0i(x,O)0ixx(x,O) - 

^ 2 (x,0)] + 0(t3) 

(x 2 (x - |) 2 < £) (2k) 



/ 
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The q uant ities 0-j_, 0^_, and their derivatives required to evaluate the 

functions F, F, and f as well as the particular integral may he 
found from the equations (15) and (l6). These quantities can he expressed 
in series form with the use of the relations (22). 


The range of integration for the integral containing f ( I ) should 
extend from -« to The series employed, however, does not converge 

over this entire range. The contribution of this integral is small 
outside the range -1 to 1; hence, with good approximation, the inte- 
gration may he taken from -1 to 1. In some cases the contribution 
of the integral involving f(l) is so small that it may he neglected 


altogether. In this case, equation (2^) converges for — lx - | I < it. 

2pn ■ 9 

Since th e maxi mum value of x ■ | is the chord, the series form for the 
correction is applicable to cases where the chord is less than /l - Mco^ 
times the t unn el height. When the integral involving f(|) is considered, 
the series solution is applicable to cases where the chord is less than 

— \Jl - Mo ,, 2 times the tunnel height. 


The constriction effects for an open tunnel may he found by an 
analysis similar to that given for the closed tunnel. Here the boundary 
condition of zero normal velocity at the wall is replaced by the con- 
dition of constant pressure along the jet boundary. To the first order 
the jet boundary is located at y* = h* . However, for the second-order 
solution the boundary condition must be satisfied on the jet boundary. 
The location of this boundary for the second-order solution may be 
determined from the first-order velocities at y* = h* . The fact that 
the ordinate of the jet boundary, is not a constant and is initially 
unknown for any Iteration higher than the first makes the solution of 
this problem more tedious than that for the closed tunnel. 


TURREL-WALL CORREGTIORS FOR A SYMMETRIC PARABOLIC-ARC AIRFOIL 


As an example of the use of the equations developed herein, the 
tunnel- wall corrections are evaluated for a parabolic-arc airfoil. This 
profile proves convenient since the first-order wall corrections can be 
obtained in terms of tabulated functions. Thus the rapidity of convergence 
of the alternating- series form of solution is readily ascertained by 
comparison with this solution. 

The parabolic-arc airfoil is defined by the equation 
Y(x) = tYq(x) = t(l - x2) 


= 0 


(-1 < x < 1) 
(x > 1; x < -1) 



NACA TN 2350 


15 


From equation ( 16 ) the u- component of velocity on the airfoil in an 
unbounded stream is given to the first order as 






x - i 


d£ 


= i ( 2 + iloge |rtfl) 


( 25 ) 


where denotes the Cauchy principal value. From equation ( 15 ) the 
u-component of velocity on the airfoil in the t unn el is 

up = ~ 1 Y l , ^) coth Mx- I) dl =- S coth X(x - l)d£ 

and, after an integration by parts, 

Uq = - J logg | sinh X(x - I) | di + 

— jlogg |sinh l(x - l) | + logg jsinh X(x + l) 

Writing the hyperbolic sine in terms of exponentials, making a change 
of variable, and performing the elementary integrations gives 

% = " P f 1 ° - Sp - |V| dv + f||^°8e I 8ial1 X ( x “ !) | + 

a w - 1 1— 

logg | sinh X(x + l) | - 2 lx + 2 logg 2 J 
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where 

v = 1 - 

a = 1 - e - a ( x+:L ) 

1, . 1 - e- 2 ^*- 1 ’ 


The Spence integral is defined by 

~>x 


m(x) , f to . 

J i X' - 1 


Thus with the use of this integral and equation ( 25 ) the tunnel- wall 
correction to the velocity is 


^i-ro5 l W 


logg | sinh X(x - l) 


+ 


log e sinh X(x + l)J - 2Xx + 2 logg 2 - 2 - x logg 


1 

1 



(2 6 ) 


The Spence integral is tabulated for positive values of x in 
reference 11. The limit b is negative for all values of- x for 
1^0. However it is easily shown that the value of the Spence integral 
for negative values of the argument is given by the relation 


RZ(-x) = ^RZ(x 2 ) - RZ (x) +|RZ(0) 

The first-order solution for the tunnel-wall correction in series 
form is found from equation ( 26 ) as • 


00 




Kx - 4) 2n_1 d| 



(-D 


n-1 


2 2a X 2a B n 
(2n)In(2n + l) 
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Similar series expressions for 0p, and- their derivatives are 

easily found. Then, -with the use of these first-order solutions, the 
functions F(t )_ and F( £) may he evaluated. The integral involving 
the function f(i) is very small and may he neglected. The tunnel- 
wall correction to the second order is then found from equation (24). 


RESULTS AND DISCUSSION 


The first-order constriction correction to the x-conqponent of 
velocity at the surface of a symmetric paraholic-arc airfoil in a wind 
tunnel is shown in figures 2 and 3 for several ratios of tunnel height 
to chord. It appears, at least for the paraholic-arc airfoil, that 
the error made in neglecting the variation of interference-velocity 
correction along the chord is small provided the chord is much less 
than the tunnel height and the Mach number is small. For these cases, 
the constant term of equation (23) gives a good value for the correction. 
However, these are just the cases for which the correction is so small 
that it is often neglected. For cases where the correction is signifi- 
cant (large Mach numbers and/or small ratios of tunnel height to chord), 
the error made in neglecting the variation of velocity correction along 
the chord may he 10 to 15 percent. Comparison of the correction (fig. 2) 
given by equation (2 6) with that given by equation (23) indicates good 
convergence of the series solution. The proper value of the correction 
is zero at a stagnation point. The fact that the correction does not 
approach zero at the stagnation points gives another 'example of the 
inability of this approximate theory to represent adequately the flow 
in these regions. 

A series for the first-order correction to the velocity is given 
by equation (23); the first term gives the often-quoted result that the 
tunnel-wall correction varies as 1/(3^ whereas succeeding terms contain 
other powers of fi. The constriction correction to the second-order for 
the parabolic-arc airfoil at midchord is shown in figure 3* It is of 
interest to note that the correction, to the second order, at the mid- 
chord varies with Mach number approximately as l/P^, the second-order 
correction being about 20 percent higher for t = 0.1 and about 40 per- 
cent higher for t =0.2. It therefore appears that the second-order 
correction should be considered for test conditions where the first- 
order correction is significant. 

The constriction corrections commonly employed to correct tunnel 
force or pressure data are applied to the stream velocity and hence to 
the dynamic pressure, Mach number, and so forth since the correction has 
a constant value over the chord. On the other hand, the corrections 
developed herein vary over the chord of the airfoil and consequently 
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■they must "be employed locally to correct the surface velocity or pressures. 
These second-order constriction corrections are more accurate them those 
normally employed hut their calculation and use entails more labor. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., February 2, 1951 
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Figure 3- - Variation of first- and second-order velocity correction at 

the midchord with Mach number. 
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